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Abstract 



We present results of molecular dynamics (MD) simulations for a non- 
entangled polymer melt confined between two completely smooth and repul- 
sive walls, interacting with inner particles via the potential f7 wa ii = (o~/z) 9 , 
where z = |-z pa rticie — z W aii| and a is (roughly) the monomer diameter. The 
influence of this confinement on the dynamic behavior of the melt is studied 
for various film thicknesses (wall-to-wall separations) D, ranging from about 
3 to about 14 times the bulk radius of gyration. A comparison of the mean- 
square displacements in the film and in the bulk shows an acceleration of 
the dynamics due to the presence of the walls. This leads to a reduction of 
the critical temperature, T c , of mode-coupling theory with decreasing film 
thickness. Analyzing the same data by the Vogel-Fulcher-Tammann equa- 
tion, we also estimate the VFT-temperature Tq(D). The ratio Tq(D) /Tq^ 
decreases for smaller D similarly to T C (D)/T c bulk . These results are in qual- 
itative agreement with that of the glass transition temperature observed in 
some experiments on supported polymer films. 
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I. INTRODUCTION 



Polymer science has had a major impact on the way we live. Just 50 years ago, materials 
we now take for granted were non-existent. Due to their structural complexity, polymers 
are generally not crystalline at low temperatures. Rather they exhibit an amorphous, glassy 
structure. The concept of the glass transition thus plays an important role in understand- 
ing the properties of polymer systems. Polymers are often used as protective coatings in 
microelectronics • In such applications, the polymer is confined in a film geometry. An 
important information for materials design is therefore how thermal properties of a polymer 
system are affected by the film geometry, in particular, whether and how the glass transition 
temperature, T g , is influenced by confinement. 

In addition to its technological importance, the investigation of the glass transition in thin 
polymer films is also from a great theoretical interest. This is closely related to the unresolved 
nature of the glass transition. Despite considerable experimental and theoretical efforts, 
there is still not a fully satisfactory description of this phenomenon. Phenomenological 
theories of the glass transition, such as the free volume theory [§-0], are attractive as they 
give a simple description of the basic observations but they contain adjustable parameters 
whose physical significance is unclear. Thermodynamic approachs, like Gibbs-DiMarzio 
theory , treat the glass transition as an ordinary thermodynamic phase transition but are 
faced with the problem that the extrapolation of the theoretical results leads to unphysical 
predictions like the existence of a temperature range for which the system entropy becomes 
negative (Kauzmann paradoxon). 

The so-called mode-coupling theory (MCT) is perhaps the most sucessfull of all descrip- 
tions of the glass transition. Within this theory, there is a critical temperature, T c , at which 
the (structural) relaxation times diverge while the static properties of the system remain 
liquid-like. The system freezes at T c keeping its amorphous structure. Thus, from the point 
of view of the MCT, the glass transition is a purely dynamic phenomenon [p!0|-[T^ . 



However, it must be stressed that neither MCT nor other descriptions of the glass tran- 
sition have a definitive character. One is therefore interested in finding ways to verify the 
basic concepts of different approachs. 

Let us assume that the slowing down of the dynamics when approaching T g can be 
interpreted in terms of the critical slowing down of the dynamics near a second order phase 
transition |15| . A second order phase transition is usually characterized by the divergence of 



a length scale. The idea of a diverging length scale related to the glass transition has given 



rise to the concept of "cooperative motion" empirically introduced by Adam and Gibbs [16 



According to Adam and Gibbs, near the glass transition, individual particle motion is frozen 
out. Thus, the only possibility for structural relaxation is that of the collective motion of 
many particles. Note that, here, the focus is no longer on the static correlations, but on 
correlations between the dynamics of particles. The associated length scale is thus a dynamic 
one. Let £(T) denote the typical size of a cluster of cooperatively moving particles. It is 
well-known that close to a second order phase transition the relaxation time of such a cluster 
scales like £ z , where z>0 is the so called dynamic critical exponent jI7||. The sharp rise of 
the relaxation times near T g is then explained by assuming the divergence of £ when lowering 
the temperature towards T g . 

The above reasoning, however, is based on an empirical assumption that £(T) increases 
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with decreasing temperature. A significant improvement was achieved by Edwards and 



Vilgis fll8fl . These authors introduced an exactly solvable model system exhibiting glassy 
behavior at low T and showed that the concept of cooperative motion alone was enough to 
give raise to a Vogel-Fulcher-Tammann (VFT)-law. 

Recent computer experiments also support the idea of cooperative motion. For example, 



it was observed by Kob et al. [|19| , pQ|| that particles move mainly in string-like clusters. 
Bennemann et al. [BTJ report on a growing length scale for the dynamics of a polymer 
melt. Strong heterogeneity in the dynamics has also been observed in molecular dynamics 
simulations of bond breakage processes and of binary mixture of soft spheres [P5|-|2"5|] . It 
is found that particles move preferably within mobile clusters thus leading to a heterogeneity 
in the dynamics. 

While easily detectable in a simulation, the regions of heterogeneous dynamics are un- 
fortunately not as easily accessible to real experiments. This results from the fact that 
dynamic heterogeneity is not necessarily strongly correlated to density fluctuations. The 
structures of these regions are therefore more or less identical. Thus, one cannot use scat- 
tering experiments to determine the length scale £(T) ||26|| . Fortunately, there is an issue 
to this problem. Recall that the relaxation time of a cluster of strongly correlated particles 
such as that observed in cooperative motion scales like £ z , where z > 0. As the temper- 
ature decreases £(T) becomes larger and eventually reaches the system size, £(T) = L. If 
this occurs, the relaxation time, r, of the system will scale like r oc L 2 . In other words, 
the relaxation dynamics of the system will become size dependent. This size dependence 
is indeed observed in Monte Carlo studies of the so called bond-fluctuation model (BFM) 
in 2D ||27|| . In the mentioned reference, the system size was varied while maintaining the 
periodic boundary conditions (pbc). An acceleration of the dynamics of the smaller systems 
have been observed in accordance with r oc L z . 

A simple way of changing the system size in a real experiment is, for example, to vary 
the thickness of a planar film. Applying the same arguments as given above to a thin film of 
thickness D we should expect finite-size effects on the dynamics for temperatures at which 
£(T) « D. Note that boundary conditions are no longer periodic, but can change from an 
absorbing one to an approximately neutral or a repulsive one. 

In this context, experiments on (model) systems reveal a stratified phenomenology. If 
the interaction between the polymers and the substrate is attractive, the glass transition 



temperature T g of the films becomes larger than the bulk value for small film thicknesses p8 



Intuitively, this effect can be attributed to chains which are close enough to the substrate 
to 'feel' the attractive interaction. The motion of these chains should be slowed down with 
respect to the bulk. In a thin film almost all chains touch the attractive substrate. So, T g 
should increase. 

On the other hand, measurements (by ellipsometry) of T g for polystyrene (PS) films (of 
rather large molecular weights) on a silicon substrate showed a significant decrease of T g from 
375K down to 345K for the smallest film thickness of lOnm, i.e. a relative change of 10% 
in T g was observed El. There have also been many experiments in recent years on freely 
standing polystyrene films (i.e., no solid substrate, but two polymer-air interfaces) [|30H 3"3H 
showing a dramatic decrease of T g by up to 20% if the film thickness becomes much smaller 
than the chain size. An interesting explanation of this observation in terms of an interplay 
between polymer-specific properties and free- volume concepts has been proposed pi. This 
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decrease becomes much weaker if one or even both of the free interfaces are replaced by 
a weakly interacting solid substrate. Whereas the strong depression of T g in the freely 
standing film could possibly be attributed to the significant release of geometric constraints 
at the air-polymer interface, the acceleration of the structural relaxation of a polymer melt 
between two (almost) neutral solid substrate is much harder to understand intuitively. 

Recent simulation results emphasize the fact that the increase or decrease of the glass 
transition temperature does strongly depend on the interaction between the polymer chains 



and the substrate |35 |. For a model of square- well spherical interaction sites interconnected 



by fully flexible strings one observes a reduction of T g in the case of a weak attraction 
between the substrate and the chains, whereas for the case of strong attraction an increase 
of T g is found fl35] . 

We study a continuum model where polymer chains are confined between two identical, 
ideally flat and purely repulsive walls [p6| , |37| , obtaining complete information in both space 
and time in atomistic detail. Instead of focusing on a computation of T g , we rather check 
some predictions of the mode-coupling theory and investigate the dependence of the mode- 
coupling critical temperature, T c , on film thickness D. We also examine the -D-dependence 
of the VFT-temperature, T , with the result that the ratios T (D)/T bulk and T C (D)/T c bulk 
behave similarly. 

After a presentation of the model in the next section, we discuss in section [III] the 
reliability of the system dynamics obtained from molecular dynamics (MD) simulations 
within NVT and NpT ensembles. In section |V| we present MD results on the influence of 
the walls on the system mobility and its critical behavior at low temperatures for a variety of 
film thicknesses showing that the critical temperature is lowered for stronger confinements 
(smaller film thicknesses). Section is devoted to a brief analysis of the local dynamics 
where we address a subtle point concerning the proper definition of local quantities in terms 
of particle positions and/or momenta. A conclusion compiles the results. 



II. MODEL 



We study a Lennard- Jones (LJ) model for a dense polymer melt |38],[39| of short chains 



(each consisting of 10 monomers) embedded between two completely smooth, impenetrable 
walls 



6,37 



Two potentials are used for the interaction between particles. The first one is 
a truncated and shifted LJ-potential which acts between all pairs of particles regardless of 
whether they are connected or not, 



U L j{r) - U hi (r c ) if r < r c , 
otherwise 



where U^j(r) =4e (o~/r) 12 — (o~/r) 6 and r c = 2 x 2 1 / f V. The connectivity between adjacent 
monomers of a chain is ensured by a FENE-potential 



k 

f^FENE(r) = ---Roln 



R 



(1) 



where k = 30e/a 2 is the strength factor and Rq = 1.5a the maximum allowed length of a 
bond. The wall potential was chosen as 
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cx\ 9 



Uw(z) = e[-) , (2) 

where z = |%> ar ticie — ^waii| (-Zwaii = ±-D/2). This corresponds to an infinitely thick wall 
made of infinitely small particles which interact with inner particles via the potential 
45e(cr/r) 12 /(47rp wall cr 3 ) where p wa n denotes the density of wall particles. The sum over the 
wall particles then yields e(o~/z) 9 . All simulation results are given in Lennard- Jones (LJ) 
units. All lengths and energies are measured respectively in units of o and e, temperature 
in units of e/k-Q (/cb = 1) and time in units of (ma 2 /e) 1 / 2 . 

The upper panel of Fig. [I] compares the bond potential, i.e. the sum of LJ- and FENE- 
potentials, with the LJ-potential. It shows that the bonded monomers prefer shorter dis- 
tances than the nonbonded ones. Thus, our model contains two intrinsic length scales. Since 
these length scales are chosen to be incompatible with a (fee or bec) crystalline structure, one 
could expect that the system does not crystallize at low temperatures, but remains amor- 
phous. This expectation is well borne out. The lower panel of Fig. [I] the structure factor, 
S(q), of the bulk system with that of a film of thickness D = 10 at T = 0.46 (note that the 
mode-coupling critical temperature of the present model in the bulk is T c bulk = 0.45 f4(J). For 
the film, two structure factors are chosen, corresponding to two different regions: a layer of 
thickness 3a centered in the middle of the film (called "film center") and the region between 
this layer and the walls (called "near walls"). In both cases, S(q) is calculated parallel to 
the walls [i.e. q — \q\, q = (q x , q y )} by averaging over all monomers in the respective region. 

Fig. [I] shows that the structure of the melt in the bulk and in the film is characteristic 
of an amorphous material. At small q, S(q) is small, reflecting the low compressibility of 
the system. Then, it increases and develops a peak at q max which corresponds to the local 
packing of monomers (27r/g max ^l) before it decreases again and begins oscillating around 
1, the large-g limit of S(q). 

The most prominant differences between the bulk and the film are found for small q and 
for g max . The compressibility of the film is higher, the value of g max is shifted to slightly 
lower q and the magnitude of S^max) is smaller than in the bulk. In the bulk one can 
observe a similar shift of g max and decrease of 5'(g max ) as the temperature increases. Also 
the compressibility of a bulk fluid increases with temperatures. Therefore, the local packing 
of the monomers in the film seems to resemble that of the bulk at a higher temperature. 
Since the local structure of the melt has an important influence on its dynamic behavior in 
the supercooled state [IT|, Fig. [1] suggests that the film releases more easily than the bulk 



at the same temperature. Indeed, we will see later that the dynamics of the system is much 
faster in the film than in the bulk when compared at the same temperature. 

All simulations have been carried out under constant normal pressure -PN.ext = p = 1- 
However, to adjust the normal pressure, we do not vary the wall-to-wall separation, D, but 
the surface area. For each temperature, the average surface area is calculated by an iterative 
approach Ell]. The system is then propagated untill the instantaneous surface area reaches 
the computed average value. At this point the surface area (and thus the volume) is fixed 
and a production run is started in iV^/T-ensemble, where the system temperature is ad- 



justed using the Nose- Hoover thermostat [f42| , f43l| . More details about the applied simulation 



techniques can be found in |^J4l|j44] . 



5 



III. NVT AND NPT VERSUS NVE: INFLUENCE ON THE DYNAMICS 



It was mentioned in section [Tl] that the production runs were performed at constant 
volume and temperature using Nose- Hoover thermostat. This thermostat slows down or 
accelerates all particles depending on the sign of the difference between the instantaneous 
kinetic energy of the system and the desired value given by the imposed temperature, i.e. 
3Nk~BT/2 (N is the number of particles) P2|]%3 ,35-fl£[|. One may therefore ask how reliable 
the resulting dynamics is when compared to pure Newtonian dynamics in the microcanonical 
(NVE) ensemble. This question was already examined for the bulk system in [39]. In this 
case, the results obtained within constant energy (NVE) simulation and that using the 
Nose-Hoover thermostat are identical. 

We are going to show that the presence of the walls does not change this behavior. 
However, we go a step further and also investigate the influence of volume fluctuations on 
the system dynamics. This point is very important if one is interested in constant pressure 
simulations. We will see that, contrary to the case of Nose- Hoover thermostat, the system 
dynamics is strongly perturbed when the system volume is allowed to fluctuate (therefore, 
for a given normal pressure, the corresponding average volume was first computed for each 
T in a ApT-simulation and then the dynamics was analyzed in production runs at constant 
volume). For this purpose, we compare results obtained from MD simulations within the 
NVE (microcanonical), NVT (canonical, using Nose- Hoover thermostat) and NpT (Nose- 
Hoover thermostat plus the fluctuations of the surface area) ensembles. 

Note that all results to be discussed in this section correspond to a film of thickness 
D = 5. Recall that D stands for the wall-to-wall separation. The distance of the closest 
approach of a monomer to a wall is approximately its own diameter (i.e. a = l). Therefore, 
a value of D = 5 corresponds to the extreme case of three monomer layers only. 

Let us first consider the velocity autocorrelation function defined as 

N N 

a(t) = (E^(°x*))/(E^(°)}- ( 3 ) 

i i=l 

Figure |2] presents results for C v (t) obtained from simulations in the NVE, NVT, and 
NpT ensembles. As seen from this figure, no difference is observed for C v (t) within various 
ensembles. The velocity autocorrelation function vanishes so rapidly that it can be equally 
well computed within all these ensembles. However, quantities that evolve slowly in time 
exhibit a different behavior. An example is the mean-square displacement (MSD) of chain's 
center of mass, 

1 - 

^3,||W = ^E<[^c m , l ,||(t)- J Rc m ,i,||] 2 > • (4) 
i=l 

Here, M is the number of chains and -R C m,i,|| is the projection of the center of mass vector of 
the i-th chain onto a plane parallel to the wall. Figure |3| shows that data for g%\\ (t) within 
NVE and NVT ensembles are identical. However, the result obtained from the NpT- 
ensemble differs strongly from the reference NVE-curve. This discrepancy is due to the 
relative small box size of a typical MD simulation. Note that the relative volume fluctuations 
scale as 1/yN. While fully negligible for real systems (N « 10 23 ), these fluctuations become 
important in a simulation where the particle number is of order of 1000. 
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IV. DEPENDENCE OF T c ON FILM THICKNESS 



We now focus on the influence of the walls on the sluggish dynamics of the system. 
For this purpose, it is instructive to recall some important features of the present model in 
the bulk at low temperatures. In Fig. 01 the mean-square displacements of the innermost 



1 M 

9l(t) = JfYl (KmncrW ~ pinner (0)] 2 ) , (5) 



1=1 



and of the chain's center of mass, gs(t), are displayed versus time for T = 0.48. For short 
times, the motion of the system can be described by assuming free particles (ballistic regime): 
gi(t) = {v 2 )t 2 = 3Tt 2 and gs(t) = 3Tt 2 /N p , where (v 2 ) is the mean-square monomer velocity 
and iVp = 10 is the number of the monomers of a chain. In agreement with the predictions 
of the mode-coupling theory (MCT) [jlOHT^I , a plateau regime emerges after the ballistic 



motion. The tagged particle "feels" the presence of its neighbors and, as the temperature 
is low, remains temporarily in the cage formed by these neighbors. However, contrary 
to a simple (atomic) liquids, where a direct crossover from the plateau into the diffusive 
regime occurs, an intermediate subdiffusive regime emerges due to the connectivity of the 
monomers ||49|| . In this regime, the center of mass already crosses over to the asymptotic 
diffusive motion, gs(t) ~ t, whereas the motion of the innermost monomer is described by a 
power law gi(t) ~ if with an effective exponent x ~ 0.63. The innermost monomers reach 
the diffusive limit only if gi(t) is larger than the end-to-end distance of a chain and thus 
outside of the time window shown in the figure. 

In Fig. [|we compare the mean-square displacement of the innermost monomer gi(t) for 
films of various thicknesses with that of the bulk. The upper panel corresponds to a high 
temperature of T = 1 where the system properties are liquid-like. At this temperature the 
influence of the walls is rather small so that gi(t) of the bulk almost overlaps with that of 
the film for thicknesses D > 10. However, the lower panel of Fig. |5] shows that the effect of 
the walls on the mobility becomes significant at all studied film thicknesses with progressive 
supercooling. Outside the initial ballistic regime, the motion is faster, the smaller the film 
thickness. 

To quantify this observation, we define relaxation times as the time needed by a given 
mean-square displacement (like gs, the MSD of chain's center of mass or g±, the MSD of the 
innermost monomer) to reach the monomer size 

gi(t = r) := 1 (defining equation for r). (6) 

Note that, due to the film geometry, there are only two independent directions (each of them 
parallel to the walls) for diffusive motion compared to three in a homogeneous melt. For 
the film, we thus compute using these two directions and then multiply the result by a 
factor of 3/2. This multiplication is necessary if the film data are to be compared with the 
corresponding bulk quantities. 

Using Eq. (H), we computed r(gi = 1) as a function of temperature for various film 
thicknesses, where, in addition to g% and gx, the MSD of all monomers, 
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1 N 

00(f) = ^£<[n(f)-n(o)] 2 ) 

i=l 



(7) 



and that of the end monomers (chain ends), 



1 M 

04(f) = ^ £ ([ r *.<=nd(f) - r ifiD d 



(8) 



have also been used. 

For each film thickness, this yields four different relaxation times as a function of temper- 
ature. As the mode-coupling theory [|ll|,|l^,|l^,|l^,0] has been rather successful in describing 
the slow dynamics of the present model in bulk HOJSOJSTf] , we tried to fit r(T) via a power 
law, 



t(T)<x\T-T c {D)\ 



- 7 (D) 



(9) 



Such a power-law divergence of the a-relaxation time is predicted by MCT for the 
bulk O,0,H. 

The fit is done as follows: First, all (three) parameters of Eq. (Q) were adjusted. The 
values of 7 obtained from r(g = 1), r{g\ = 1), r(g 3 = 1) and r(g 4 = 1) agreed well within 
the error bars. Therefore, we fixed 7 at the average value for the given film thickness and 
repeated the fits. 

Table | contains results for T C (D) and j(D) obtained in this way. Figure |B] shows a 
representative example for this analysis. It depicts r -1 / 7 versus T for a film of thickness 
D = 5 (note that D is the wall-to-wall distance, the thickness of the region with non- 
vanishing monomer density is approximately D — 2, see Fig. [19]). The intersection of the 
straight lines (MCT-fit results) with the T-axis yields the critical temperature at this film 
thickness: T C (D = 5) = 0.305 ± 0.005. Note that, despite the highly non-linear relationship 
between the MSD's used to define the various r's, all fits yield the same T c . 

To test this analysis the resulting critical temperature can be used to linearize reduced 
plots of the relaxation times versus T — T c on a log-log scale. Figures [7| and |8| show that 
the power law (|9|) motivated by the ideal MCT is a good approximation of the data at 
temperatures close (but not too close) to T c (see, for example, Fig. [K] and also |40| , |50|j5^ 
for comparable bulk data). 

However, as indicated by the solid line in the both figures, r(T) can also be described 
by a Vogel-Fulcher-Tammann (VFT)-fit in the studied temperature range, i.e. by 



t(T) oc exp 



c{D) 



T-T (D) 



(10) 



where c is a constant which can depend on film thickness. The possibility of describing the 
same data both by a power law (MCT) and by a VFT-fit has also been observed for bulk 
properties of the present model (see Fig. 10 in |J?J). We therefore use the VFT-formula as an 
independent approach and determine the VFT-temperature, T , for various film thicknesses. 
Table I contains the results for T (D) thus obtained. A plot of T C (L>) /T c bulk and r (D)/T bulk 
is shown in Fig. R. Both the mode-coupling critical temperature and the VFT-temperature 



S 



exhibit similar D-dependences, thus suggesting that the presence of the smooth walls results 
in a reduction of the glass transition temperature, T g . As already mentioned in section |, a 
reduction of the glass transition temperatures has also been reported both from experiments 
on supported polymer films |29| and on free standing polymer films and also from 



MD simulations of a model of square-well spherical interaction sites interconnected by fully 
flexible strings in the case of a weak attraction between the substrate and the chains ||35|| . 
In Fig. ITUl r(<7i = 1) is depicted versus T—T c for the homogeneous melt (bulk) and films 



of various thicknesses. The solid line indicates the power law r ~ (T — T c bulk ) 7bulk with 
T c bulk = 0.45 and 7b u ik = 2.09 P0| , |5D| , |5T| . To show the (slight) .D-dependence of the exponent 



7, a fit to D — 5 data is also depicted (see the long-dashed line with an slope of 2.5). It is seen 
from this figure that both in the film and in the bulk, for temperatures very close to T c , the 
relaxation times increase more slowly than predicted by the ideal MCT. This discrepancy is 
an indication of slow relaxation processes which are not taken into account within the ideal 
MCT. As temperature decreases, the contribution of these processes becomes important and 
the ideal MCT no longer holds. A common picture in describing the slow dynamics of the 
system at these temperatures is that of an energy landscape with finite potential barriers. 
Relaxation processes are then described as jumps (hoppings) between neighboring energy 
minima. If one assumes a sharp distribution of the potential barriers around a given value, 
Eq, i-e. if all such barriers are assumed to have more or less the same height, the probability 
of a jump over a barrier is then equal to the probability of having an energy E > Eq. 
Assuming that the diffusion constant is proportional to this probability, one obtains, 

Diffusion Constant oc / exp(— ) dE = exp(— ). (11) 

Je k B T k B T 

Thus, due to the above picture, at low enough temperatures the diffusion constant obeys an 
Arrhenius law. Here, "low enough" means temperatures very close to or below T c . However, 
for higher (but not too high) temperatures one again expects the validity of a power law 
also for the diffusion constant, 

Diffusion Constant oc |T - T C (L>)| 7(D) . (12) 

Using the mean square displacements, we compute the diffusion constant from 

Diffusion Constant = lim — - . (13) 

Note that, in the diffusive limit (t — > oo), there is no difference between g , gi , g 3 or g 4 
(see for example the upper panel in Fig. [L4], where g (t) and g 3 (t) coinside for large t). In 
the praxis, however, one has to evaluate Eq. ( [T~3[ ) at large but finite t. As the mean-square 
displacements of chain's center of mass reaches the diffusive limit faster than other MSD's 
(see Fig. £|), one should use g 3 (t) in an evaluation of Eq. (|i~3|). 

best results on the diffusion constant one obtains it is numerically A log-log plot of the 



diffusion constant versus T— T c is depicted in Fig. [IT] for various film thicknesses and for the 
bulk. Similar to the behavior of the relaxation time discussed above, at temperatures very 
close to T c , results for the diffusion constant deviate from the power law given by Eq. (|i~2|). 



Motivated by the discussion which led to Eq. (ITT]), we try to fit the diffusion constant at 
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law temperatures by an Arrhenius law. This is demonstrated in Fig. [12] for the case of 
a film of thickness D = 10, where Eq. (|TT|) is applied to the last few data points in low 
temperature regime. From this fit, we obtain E pa 9.35. However, one should be carefull 
with an interpretation of Eq. Recall that Eq. ([11]) was derived assuming that all potential 
barriers have the same height of Eq. Therefore, an interpretation of Eq as the typical height 
of the potential barriers in an energy landscape requires a study of the inherent structure of 
the system. 



It is also shown in Fig. [L2| that the data at higher but not too high temperatures are 
well described by a power law, as expected from the MCT analysis of the relaxation times 
presented in this section. 

Next we look at the temperature dependence of a crossover time, r co , given by 

gi{t = r co ) := 0r 2 sc = 0.054 (definition of r co ) . (14) 

This definition is motivated by an analysis of the bulk, where r sc (pa 0.1 monomer diameter) 
appears as a relevant length scale for the dynamics of a monomer in its local environment 
'cage") PD||5D[ . The value 6rf c (roughly) coincides with the inflection point of g\{t) and 



thus marks the crossover from the early time regime, where gi(t) is almost flat, to the late 
time regime, where g\(t) gradually increases with time and finally becomes diffusive (see 



Fig. [L4|). Motivated by ideal MCT, one expects that r co diverges at T c with a power law, 



-1 

r co (T)cx|T-T c |^, (15) 

with < a < 0.5. The exponent a describes the initial decay of the correlation functions 
from the ballistic regime to the plateau [fL^-141- The analysis of the bulk data showd that 



a = 0.352 g. 

Figure [13] displays r c ~ 2a versus T—T c for the bulk and for various film thicknesses: D = 5 
(pa3.5R g ; i2 g = bulk radius of gyration), D = 7 (pa5R g ), D = 10 (pa7R g ) andL> = 20 (paUR g ). 
The critical temperatures are taken from Table | and for the exponent a we used the bulk 



value a = 0.352 |40|| . It is seen from Fig. [13] that this value yields a good linearization of 



the bulk data. For r co obtained from the film data, however, deviations occur indicating a 
(rather slight) D-dependence of the exponent a. 

The previous analysis of the T-dependence of the relaxation time r and of the crossover 
time r co suggests that the critical temperature depends strongly on film thickness, but 
the variation of the exponents a or 7 is rather weak. Within the framework of MCT the 
exponents a and 7 are related to one another and determined by the so-called exponent 
parameter A, which controls the time evolution of the monomer mean-square displacement 
in the plateau regime (see Fig. |J). Thus, the weak variation of the exponents a and 7 with D 
should imply that the time dependence of the mean-square displacements in the film closely 
follows that of the bulk, in the plateau regime. It is therefore interesting to test whether the 
mean-square displacements of the film and of the bulk obey the same master curve when 
compared for the same reduced temperature, T—T c . This idea is first examined in the upper 
panel of Fig. [14] where the MSD's of the innermost monomers, gi(t), and of the chain's 
center of mass, gs(t), are compared at T — T c = 0.07 for the bulk and for a film of thickness 
D = 10 (as T c bulk = 0.45 and T C (D = 10) = 0.39, the bulk data at T = 0.52 are compared to 
the corresponding film data at T = 0.46). There is a good agreement between film and bulk 
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data over a long period of time, in particular for intermediate times. The lower panel of the 
figure extends this comparison to lower temperatures. Whereas the film data still closely 
agree with each other for all times simulated, this is no longer the case for the bulk data. By 
simply shifting the T-axis it is only possible to rescale bulk and film data onto a common 
master curve in the intermediate time window of the plateau. At later times, however, gi(t) 
of the bulk increases faster with t than gi(t) of the films. This difference is more pronounced 
the closer T to T C (D). 

A preliminary analysis of the incoherent scattering function also reveals the property 
that for intermediate times, the data corresponding to different film thicknesses follow the 
same master curve when compared for the same reduced temperature T — T C (D). A more 
detailed investigation and comparison with the behavior of the mean-square displacements 
is under way. 

V. LOCAL DYNAMICS 

The discussion of the previous section illustrated that the dynamics of the film is accel- 
erated with respect to the bulk. Obviously, this is a consequence of the walls. A natural 
question is therefore how the motion of the monomers depends on their distance from the 
wall. To obtain a better insight into this dependence we study the local displacements, 
i.e. the mean-square displacements measured within layers of a small thickness ("bin") at 
various distances from the wall. However, when trying to analyse local displacements the 
problem of how to assign particles to layers has to be addressed. A possible definition is to 
associate a monomer labeled i to that bin to which its initial ^-position, Zi(to), corresponds, 

N N 

g {t;z) = (^6(z-z i (t ))(r i (t)-r i (t )) 2 /^6(z-z i (t ))} . (16) 

8=1 1 

The drawback of this definition is, however, that some particles may leave this layer at 
a later time. If one adds their contribution to the MSD to that of those particles which 
allways remain in the initial region, a sort of averaging over other regions is introduced. 
Obviously, the larger the time difference t — to, the greater the probability that some tagged 
particle has left its initial layer. As a consequence, the local character of the obtained 
information becomes questionable for large times. An estimate of the time beyond which this 
averaging over neighboring regions becomes apreciable is the time at which the mean-square 
displacements in transversel direction become comparable to the half of the bin size. For 
larger times local information is gradually lost and the MSD's of different regions converge 
towards the average MSD of the whole system. This point is nicely demonstrated in Fig. |T5| 
where the MSD of all monomers in the direction parallel to the walls, go,\\{t), is shown for 
two regions of different mobility in a film of thickness D — 10 at T — 1 (high temperature 
liquid state). As <7o,]| > 1 (half of the bin size), the MSD of both regions converge towards 
that of the whole system for large times. 

To avoid this problem, a more stringent criterion must be used to compute local MSD's 
from particle displacements: The contribution of a tagged particle to the MSD of a given 
layer at a time t should be taken into account if and only if the tagged particle has been in 
the same layer for all times to < t' < t, 
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N t N t 

g (t; z) = (52I[6(z- z t (t>)) (r*(t) - n(t ))7 E II ^ " ■ W 

i=l t'=to 1 t'=t 

This correct definition is, however, computationally mor demanding. The reason lies in the 
fact that the contibution of particles which have left their initial layer at some later time, has 
to be ignored for all subsequent times. As the number of such particles increases with time, 
the statistical accuracy decreases at late times. To improve it, one needs more independent 
samples of the same system. 



In Fig. 16 we focus on another aspect of the same problem. Figure O suggests that a 



good estimate of the time at which the difference between Eqs. ( jig ) and fll7|) is no longer 
negligible is gi \\(t) = bin width/2. Here, we stress that the inverse is not necessarily true, 
i.e. the error is not necessarily negligible for shorter times. Already in Fig. |15| one can 
observe that, in the case of the layer close to the wall, deviations from the correct curve 
occcur at shorter times than the intersection point between go,±{t) and the horizontal line 
indicating half of the layer thickness. This point is more clearly demonstrated in Fig. WE 
where the MSD's of a film of thickness D = 20 at a low temperature of T = 0.46 are used 
[T C (D = 20) = 0.415]. Here, for the layer in the inner part of the system, the deviations 
are apreciable already at times which are about two orders of magnitude smaller than the 
"intersection time" . One should therefore be aware of this discrepancy when analyzing local 
quantities. 

Figure [17| depicts the local MSD of the innermost monomer, gi(t; z), for a film of thickness 
D = 20 at a low temperature of T = 0.46. The displacements are calculated parallel to the 
walls using Eq. (|i~7|) . We see from this figure that the mobility in film center is practically 
equal to that of the bulk, whereas monomers in the proximity of the wall are much mobile. 
There is a gradual transition from the two step relaxation characteristic of the bulk at this 
temperature to a smooth crossover from microscopic to free diffusive motion as the layers 
are closer to the wall. Whereas the bulk-like dynamics for the film center is plausible -the 
monomer density profile is flat there and equal to pbuik (see Fig. the continuous speeding- 
up of the dynamics in spite of the pronounced oscillations of the monomer density is less 
intuitive. It could be related to the following points: An important factor is certainly the 
wall potential. This potential is softer [see Eq. (|2])] than the LJ-potential which means that 
a monomer close to wall can further penetrate the wall than a monomer can approach the 
center of mass of its nearest neighbors. This should yield a higher mobility of the monomers 
at the wall at low T where the dense local packing of the monomers is responsible for the 
mutual blocking and slowing down of the bulk dynamics. 

On the other hand, the static structure factor of Fig. |l] shows that the first maximum is 
not as large in the film as in the bulk at the same low temperature. This implies that the 
overall arrangement of the monomers in the nearest neighbor shells is not as pronounced in 
direction parallel to the wall as it is in the bulk. Since this local packing is an important 
factor for the slowing down of the structural relaxation ( Hansen- Ver let freezing criterion), 
the dynamics in the film should be faster. 

Using a local version of Eq. (^), we define a z-dependent relaxation time, t(z), 

g%(t=T(z); z) := 1 (defining equation for t(z)), (18) 

where i — 0, 1, 3, 4 denotes various types of displacements of a polymer system [see sections [TO 
andffVt 
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For a film of thickness D = 20, we applied Eq. ( pT8|) to the local MSD of all monomers, 
g Q (t; z), and computed t(z) at various temperatures. The results thus obtained are displayed 
in Fig [18]. Not unexpectedly, the effect of the walls is quite small at high temperatures. At 
T = l, for example, there is a wide region of z- independent relaxation time around the film 
center. This constant value agrees with the relaxation time r obtained by applying Eq. (||) 
to the bulk data. As the temperature decreases the presence of the walls is "felt" also in 
the inner part of the film and the width of the region of constant relaxation time decreases. 

This propagation of the wall effects into the inner part of the film is also observed in 
the density profile. Figure ITU] shows that density oscillations become more pronounced and 
long ranged at lower temperatures. Since the normal pressure is kept constant during the 
simulations (as it was also the case in preceeding simulations in the bulk ^(|), the average 
value of the density in the film center increases with decreasing temperature and coincides 
with that of the bulk at the same temperature. This increase of the density in the film gives 
raise to more pronounced oscillations of the monomer profile. 



VI. CONCLUSION 

We presented results on molecular dynamics simulations for a model of non-entangled 
short polymer chains confined between two perfectly flat, non-adsorbing and impenetrable 
walls. The monomer-wall interaction is modeled by potential which diverges as z~ 9 when 
a monomer approaches the wall. This repulsion is weaker than the 12-6 Lennard- Jones 
potential of the monomer-monomer interaction and also than the soft-sphere potential (~ 
z~ 12 ) used in simulations of confined supercooled simple liquids This special choice 



of the monomer-wall interaction has an influence on both static and dynamic properties of 
the polymer films. 

We find oscillations of the monomer density which start from a large value close to the 
wall and decay towards the bulk density in the middle of the film for sufficiently large film 
thickness. The amplitude of these oscillations becomes more pronounced with decreasing 
temperature, since the bulk density increases in our constant-pressure simulations. How- 
ever, the height of the largest maximum at the wall does not exceed 1.5 times the bulk 
density in the temperature range studied. This is much smaller than the difference found in 
simulations where the particle-wall interaction is modeled by a soft-sphere potential. There, 
the density can be more than 4 times larger than the bulk density at low T, which leads 
to a vanishing of the density in the subsequent minimum. In this case it is conceivable 
that the smoothness of the walls, which could promote fast particle motion ( "slip boundary 
condition"), is completely outweighed by the high density in the first layer, and one can 
speculate that the subsequent layers are also slowed down with respect to the bulk. This 
could explain why the studies [p4,55| found an increase of the glass transition temperature 



rather than a decrease although the confinement was realized by completely smooth walls 
as in our work. 

In our work the repulsive monomer-wall interaction is much softer, even than that be- 
tween the monomers in the bulk. Therefore, a monomer can come closer to the wall than to 
anyone of its nearest-neighbors. This should give a monomer at the wall more freedom to 
move. Since the slowing down of the structural relaxation of our model in the bulk mainly 
results from the blocking of a particle by its nearest neighbors PD|J5^|,|57[] , one can expect 
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this difference in repulsive interaction to become particularly important at low temperatures 
close to T c . This expectation is borne out by the simulation data. Whereas the monomers 
in the film center exhibit very sluggish motion in an intermediate time interval, a signa- 
ture of the "cage effect" , this intermittence of the displacements is not at all visible for the 
monomers in contact with the wall. They behave as if they were at a higher temperature (see 
Fig. [17]). Their higher mobility also triggers faster motion of adjacent monomers, leading to 
an overall acceleration of the dynamics in the film. 

We have tried to quantify this acceleration by fitting various relaxation times and the dif- 
fusion coefficient of a chain to both a power law, motivated by mode-coupling theory (MCT), 
and by the Vogel-Fulcher-Tammann (VFT) equation. The power law yields an estimate of 
the critical temperature, T c , and the VFT-equation of the Vogel-Fulcher temperature, T , as 
a function of film thickness D. In accord with the qualitative observations described above 
both T C (D) and T0(D) decrease with decreasing film thickness. Furthermore, a comparison 
of the rescaled quantities, i.e., of T c (Z))/T c bulk and T Q (D)/TQ Ulk , shows that the variation 
with D of T C (D) agrees with that of T (D). Although the extrapolation to T is less reliable 
than that to T c , since it has to cover a larger T-range, this result suggests that the glass 
transition temperature of our model should follow the same behavior. 

The reduction of the critical temperature seems to be the dominant influence of the walls 
on the dynamics in the time window of the MCT-/3 process. For these intermediate times, 
the mean-square displacements of the films of various thicknesses and of the bulk exhibit the 
same time dependence when compared for a fixed reduced temperature T—T c . This shifting 
property of the temperature axis suggests that similar processes cause the intermittence of 
monomer motion on these intermediate time scales in the bulk and in the films. However, 
deviations between film and bulk dynamics are found as soon as the monomers escape from 
their local environment and the a-relaxation sets in. In order to understand these differences 
better we want to investigate the incoherent intermediate scattering function and related 
quantities which provide information on the dynamics on different length scales. Work in 
this direction is under way. 
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TABLES 



D 


5 


7 


10 


15 


20 


bulk 


To 


0.204 ± 0.007 


0.253 ± 0.013 


0.288 ± 0.006 


0.297 ± 0.007 


0.308 ± 0.004 


0.328± 0.008 


T c 


0.305 ± 0.006 


0.365 ± 0.007 


0.390 ± 0.005 


0.405 ± 0.008 


0.415 ± 0.005 


0.450 ± 0.005 


7 


2.5 ± 0.2 


2.4 ± 0.2 


2.1 ± 0.1 


2.2 ± 0.1 


2.1 ± 0.1 


2.09 



TABLE I. Survey of the VFT-temperature, Tq, mode-coupling critical temperature, T c , and 
the critical exponent, 7, for different film thicknesses D and for the bulk. To was determined via 
fits to Eq. ( |l~0|) both for the film and for the bulk. As to T c , we determined T C {D) from fits to 
Eq. @. T^ nlk was known from previous analyses |4C,5C,51|. This result for T c bulk is also obtained 
applying Eq. (S) to the bulk data. 
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FIG. 1. Upper panel: Illustration of the potentials of the model. The bond-potential results 
from a superposition of the Lennard- Jones (LJ) and the FENE potentials. The minimum position of 
the bond-potential is smaller than that of the LJ-potential. This incompatibility favors amorphous 
structure at low T which is confirmed by the behavior of S(q). Lower panel: Comparison of the 
static structure factor S(q) of the melt in the bulk and in the film (D = 10) at T = 0.46 (critical 
temperature of mode-coupling theory in the bulk: T c ~ 0.45 J40[|). 'Near walls' and 'film center' 
mean averages over the regions close to the walls (0 < z < 3.5) and the inner portion of the film 
(3.5 < z < 5). z is the distance of a particle from the (left) wall. The vertical dashed line indicates 
the q- value corresponding to the bulk radius of gyration {El ~ 2.09). 
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FIG. 2. Velocity autocorrelation function (ACF) as computed from simulation runs within 
the NVE, NVT and iVpT-ensembles. Neither the coupling to the heat bath nor the fluctuations 
of the volume seem to affect the behavior of this quantity. The inset shows a magnification of the 
initial behavior of the velocity- ACF (D = 5, P N = 1, T = 0.55, iV = 500). 
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FIG. 3. The mean-square displacement (MSD) of chain's center of mass in direc- 
tion parallel to the walls, ^ (t), obtained from NVE, NVT and NpT simulation runs 
(D = 5, Pjv = 1, T = 0.55, iV = 500). Obviously, the iVT/T-result is identical to that obtained within 
the iWE'-ensemble simulation. Contrary to that, the time evolution of g 3 ^ (t) in the A^pT-ensemble 
simulation is unphysical for t > 1. 
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FIG. 4. Log-log plot of the mean-square displacements of the innermost monomer, gi(t), and 
of the chain's center of mass, #3 (t), versus time for T = 0.48. The initial ballistic behaviors for gi(t) 
and 53 (t), i.e., gi(t) = (v 2 )t 2 = 3Tt 2 and 53 (i) = 3Tt 2 /N p ((v 2 ): mean-square monomer velocity, 
N p = 10: chain length), and the late time diffusive behavior are indicated as dashed lines. In 
addition, a power law fit gi(t) ~ t x with an effective exponent x — 0.63 is shown as another 
dashed line. The dashed horizontal lines represent the radius of gyration R 2 (~ 2.09; upper line) 
and the plateau value 6r 2 c of a MCT-analysis (~ 0.054; lower line), respectively. The dashed 
vertical line indicates the time, r co , of the intersection between gi(t) and 6r 2 c . Additionally, the 
mode-coupling approximation in the /3-relaxation regime is shown as a thick solid line. Adapted 
from reference [f)2^ . 
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FIG. 5. Log-log plot of the mean-square displacement of the innermost monomer gi(t) for two 
different temperatures: T = 1 (high temperature, normal liquid state; upper panel) and T = 0.5 
(supercooled state close to T c bulk = 0.45; lower panel). The figures compare bulk data with the 
displacements measured parallel to the walls in films of different thickness: D = 5 (~ 3.5i? g ), 
D = 10 (« 7.Rg) and D = 20 (« 14i? g ). The film data were multiplied by 3/2 to account for the 
different number of spatial directions used to calculate g\ (t) (i.e., 3 directions for the bulk, but only 
2 for the films) . The bulk end-to-end distance R 2 and radius of gyration R 2 , are indicated as dashed 
horizontal lines. Furthermore, the lower panel also shows the plateau value Qr 2 c of a MCT-analysis 
(~ 0.054; lowest dashed line). The solid lines represent the behavior of g\{t) expected in different 
time regimes: ballistic at short times (gi(t) ~ t 2 ), diffusive at late times (g\(t) ~ t), and dominated 
by chain connectivity for times where gi(t) > 1 = monomer diameter (gi(t) ~ t x ; x = 0.63 = effective 
exponent). 
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FIG. 6. Plot of r^ 1 / 7 versus T. The relaxation time r was determined by Eq. (^) using the 
mean-square displacements of inner, end and all monomers and of the chain's center of mass. The 
mode-coupling exponent 7 was first determined from fits to Eq. @ where all fit parameters were 
first treated as independent. As all values obtained for 7 agreed within the error bars, 7 = 2.5 ±0.2, 
we repeated the fits with 7 = 2.5. The MCT fits to the data are represented by the straight 
lines. From the intersection of these lines with the T-axis the critical temperature is determined 
T C (D = 5) = 0.305 ±0.005. 
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FIG. 7. Relaxation time r{gi = 1) for an extremely thin film of three monomer layers only 
(wall-to-wall distance D = 5). Different mean-square displacements are used for the analysis: go 
= MSD of all monomers, g\ = MSD of the innermost monomer, g% = MSD of the chain's center 
of mass and 54 = MSD of the end monomers. The long-dashed line indicates the fit using Eq. (|9|) 
motivated by the ideal mode-coupling theory. However, the solid line which corresponds to a fit 
using the VFT-law describes equally well the data [see Eq. (|l0|)]. Both fits shown here were done 
for gi(t). 
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T-Tc 

FIG. 8. Relaxation time r(gi = 1) for a film of thickness D = 10. Different mean-square 
displacements are used for the analysis: go = MSD of all monomers, g\ = MSD of the innermost 
monomer, g% = MSD of the chain's center of mass and 54 = MSD of the end monomers. The 
long-dashed line indicates the fit using Eq. (||) motivated by the ideal mode-coupling theory. The 
solid line corresponds to a fit using the VFT-law [see Eq. fllpp]. Both fits shown here were done 
for gi(t). 
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FIG. 9. Ratios T C (D) /T^ ulk and Tq(D) /Tq u1]s - versus film thickness D. The critical tempera- 
tures, T C (D), of the films were obtained from fits to Eq. (|9|). Similarly, the Vogel-Fulcher-Tammann 
temperatures, Tq(D), and Tq u ^ are results of fits to Eq. (|io|). The error bars of Tq{D) are larger 
than those of T C (D), since Tq(D) ^ T C {D) so that the difference between the lowest simulated 
temperature and the extrapolated reslut is much larger for the VFT-temperature than for the 
critical temperature. Therefore, the results for T C (D) are more reliable. 
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FIG. 10. Relaxation time r determined from gi(t) by requiring g\ (r) = 1(= monomer diameter). 
The solid line indicates the power law r ~ (T - T c bulk )-^"iK suggested by MCT for the bulk at 
temperatures above T c bulk (T c bulk = 0.45 and 7t, u ik = 2.09 @). The long-dashed line represents a fit 
to .D = 5 data with an exponent of ^{D = 5) = 2.5. Both in the film and in the bulk, for temperatures 
very close to T c , the relaxation times increase more slowly than predicted by the ideal MCT. 
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FIG. 11. Diffusion constant of a chain versus T — T c for films of thickness D = 5, 10, 20 and 
for the bulk. The critical temperatures, T C (D), of the films were obtained from fits to Eq. @ 
sec table |. T c bulk was known from the previous bulk analysis pOfl . The straight lines indicate 
fits using Eq. (|i~2]). The mode-coupling exponent 7 is larger for stronger confinement (smaller film 
thickness). Note that, the lowest simulated temperature is T = 0.35 for D = 5 and T = 0.4 for 
D = 10. However, as the critical temperature of the thinner film is much lower (T C (D = 5) = 0.305 
compared to T C {D = 10) = 0.39) the D = 10 data are closer to the corresponding T c . This explains 
why practically no deviation from a power low is observed for D = 5. 
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FIG. 12. Diffusion constant of a chain versus inverse temperature for films of thickness 
D = 5, 10, 20 and for the bulk. The critical temperatures, T C (D), of the films were obtained 
from fits to Eq. @ [see table ||. was known from the previous bulk analysis jttj. The 

solid line indicates the fit results using Eq. ( |l2|) for D = 5. The dashed line indicates a similar 
fit result for D = 10. A short solid line indicates the result of an Arrhenius fit [see Eq. (11)] to 
the D = 10 data in a low temperature range, where the ideal MCT no longer holds [see also the 
text]. Note that, the lowest simulated temperature is T = 0.35 for D = 5 and T = 0.4 for D = 10. 
However, as the critical temperature of the thinner film is much lower (T C (D = 5) = 0.305 compared 
to T C (D = 10) = 0.39), the D = 5 data are farther from the corresponding T c . This explains why a 
power low fits well to D = 5 data at all temperatures. 
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FIG. 13. Temperature dependence of the crossover time r co and the relaxation time r{g\ = 1) 
(lower panel) for the bulk and different film thicknesses: D = 5 (~ 3.5i? g ; R g = bulk radius of 
gyration), D = 7 5i? g ), D = 10 7R g ) and D = 20 (ss 14i? g ). The critical temperatures were 
obtained from fits to Eq. (|9|) [see Table |]. The crossover time was determined from the mean-square 
displacement of the innermost monomer, gi(t), by the condition <?i(t co ) = 6rg C = 0.054. This choice 



as well as power laws of both panels are motivated by a MCT-analysis of the bulk data |4C.5G,51|. 
The value 6r^ c (roughly) coincides with the inflection point of g\ (t) and thus marks the crossover 
from the early time regime, where gi(t) is almost flat, to the late time regime, where g\ (t) increases 



faster with t and finally becomes diffusive [see Fig. 14 1. The exponent a used to linearize the data 
close to T c is that of the bulk (a = 0.352 PET 
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FIG. 14. The upper panel shows the mean-square displacements of the innermost monomer, 
gi(t), and of the chain's center of mass, gs(t), between the bulk (symbols o and *) and a film 
(solid lines) of thickness D = 10 (~ 7i? g ). The panel compares the results for the same distance 
to the critical temperature T c (i.e., T — T c = 0.07). The displacements of the film are calculated 
parallel to the wall. They were multiplied by 3/2 to put them on the scale of the bulk data. The 
two horizontal dashed lines show the bulk end-to-end distance R 2 and the radius of gyration R 2 , 
respectively. The ballistic short-time behavior g\{t) = 3Tt 2 is indicated for the bulk at T = 0.52. 
The lower panel shows behavior of gi(t) for D = 5, 10, 20 and for the bulk when approaching T c 
(T C (L> = 5) = 0.305, T c (£> = 10) = 0.39, T C (.D = 20) = 0.415, T c bulk = 0.45). The dashed horizontal line 
shows the plateau value 6r^ c (~ 0.054) of a MCT-analysis performed on the bulk data [40|. For 
intermediate times, the film and bulk data corresponding to a given T — T c follow the same master 



curve. 
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FIG. 15. Illustration of the dependence of the local MSD upon the way the particle displace- 
ments are attributed to layers ("localized"). For this purpose, two regions of different mobility 
are considered: A region around the film center, 3 < z < 5, and the first layer close to the walls, 
1 < z < 3 (z — | ^particle — ^waii I ) • Note that there is practically no particle in the very proximity of 
the walls, i.e. in the region with < z < 1. If the displacement of a particle is airways attributed to 
the layer, where the particle was observed at to, the obtained results for the local MSD's are biased 
in the following sense: If a tagged particle leaves its initial layer at a later time, its motion will 
partly represent the properties of the layers it has "visited" so far. Attributing the mobility of such 
a particle to the initial layer thus represents an averaging over many layers. The error introduced 
in this way is negligible for short times. However, as the MSD in transversel direction, go,±(t), 
becomes comparable to the half of the layer thickness [see the intersection of the horizontal line 
with go,±(t)] the error dominates and the local character of the data is lost to a large extent. For 
very long times, no local information "persists" this averaging and curves belonging to different 
regions converge towards the system average [see long time behavior of the dotted lines]. The 
correct definition requires that only those particle are allowed to contribute to the mobility of a 
given region, which stay in the same layer for all times to < t' < t [sec long-dashed lines]. However, 
there is also a drawback of this correct definition: As any particle which leaves its initial layer 
must be excluded by the computation of the MSD's, the statistical uncertanity increases with the 
number of such particles and thus with time [see long-dashed lines for large t] . Statistical accuracy 
at large times thus requires more independent samples of the same system. The difference in the 
magnitude of go t ±(t) and go^^t) arises from the fact that there is only one independent direction 
perpendicular to the walls compared to two independent parallel directions. 
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FIG. 16. Same as in Fig. 15 but for a film of thickness D = 20 at a low temperature of T = 0.46 



[T C (D = 20) = 0.415]. The figure illustrates that attributing the displacement of a particle to the 
layer, where it was at the initial time to, can lead to errors which strongly depend on the layer 
under consideration. Therefore, the results obtained using this definition are not reliable even at 
times for which the MSD in transverse direction, go ±, is much smaller than the half of the layer 
thickness which, in turn, is the localization resolution. The difference in the magnitude of go t j_(t) 
and 3o,[| (*) arises from the fact that there is only one independent direction perpendicular to the 
walls compared to two independent parallel directions. 
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FIG. 17. Layer-resolved mean-square displacement of the innermost monomer, g%(z, t), ver- 
sus time for a film of thickness D = 20 (~ 14i? g ; R g = bulk radius of gyration) at T = 0.46 
(T C (D = 20) = 0.415; see Table |). z denotes the distance from the (left) wall. The displacements 
are calculated parallel to the wall and multiplied by 3/2 (to put them on the same scale as the 
bulk data (•) which are averaged over 3 spatial dimensions instead of over only 2 for the film). The 
gi (z, £)-data (indicated by lines) are averages over all monomers between some z m \ n < z < z max , 
which remain in the specified z-interval for all times shown in the figure. Due to the loss of the 
statistical accuracy at late times, the data is some times cut off at times with large statistical 
noise. In the middle of the film (9 < z < 10) g\(z,t) coincides with gi(t) of the bulk, whereas it 
is much faster at the wall (0 < z < 2). The average behavior of the film (i.e., the average over all 
layers) is shown by □. The dashed horizontal line indicates the plateau value 6r^ c (~ 0.054) of a 
MCT-analysis performed in the bulk JlC . 50|. 51 ] . 
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FIG. 18. Relaxation time t(<?q = l',z), versus distance from the wall for a film of thickness 
D = 20 (« 14i? g ; i? g = bulk radius of gyration). Here, z stands for the distance of the center of 
a layer (bin) from the (left) wall. The thickness of a bin is Az = 1. In the calculation of the 
local mean-square displacements of all monomers, go(t,z), which underlies the definition of r, we 
consider the contribution of only those monomers, which remain within the specified layer for all 
times t' <t (see also the text). At high temperatures, there is a wide region around the film center, 
where r is independent of z. Contrary to that, at low T the presence of the walls is "felt" also close 
to the film center [note that the lowest temperature is close to T C (D = 20) = 0.415; see Table [|. 
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FIG. 19. Temperature dependence of the monomer density profile p(z) for D = 20 ~ 14i? g 
(Rg = bulk radius of gyration), z denotes the position of a monomer from the (left) wall. Since 
the profiles are symmetric with respect to the middle of the film, only one half is shown. The 
temperatures are characterisitc of the high-temperature state (T = 1) and the supercooled state 
(T = 0.46 > T c (20) = 0.415; see Table |) of the melt. The horizontal filled circles (•) indicate the 
bulk densities /9buik = 0.91 and pt>uik = l-038 at T = l and T = 0.46, respectively. 
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